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ABSTRACT 

The abundance of galaxy clusters as a function of mass is determined using the 
2dFGRS Percolation-Inferred Galaxy Group (2PIGG) catalogue. This is used to es- 
timate the amplitude of the matter fluctuation spectrum, parametrised by the linear 
theory rms density fluctuations in spheres of 8/i _1 Mpc, as- The best-fitting value for 
this parameter is highly correlated with the mean matter density in the Universe, 

Cl m , and is found to satisfy a$ — 0.25 ^m°' 92 4 ' 5 '° m °' 22 ^ ± 10% (statistical) ±20% 
(systematic) for 0.18 < £l m < 0.50, assuming that Q m + Q\ = 1. This gives a$ — 0.89 
when f2 m = 0.25. A ~ 20 per cent correction has been applied to undo the system- 
atic bias inherent in the measurement procedure. Mock catalogues, constructed from 
large cosmological N-body simulations, are used to help understand and model these 
systematic errors. The abundance of galaxy groups as a function of group bj band lu- 
minosity is also determined. This is used in conjunction with the halo mass function, 
determined from simulations, to infer the variation of halo mass-to- light ratio over four 
orders of magnitude in halo mass. The mass-to-light ratio shows a minimum value of 
100]iMq/Lq in the 6j band at a total group luminosity of L&, ~5x 10 9 h~ 2 L Q . To- 
gether with the observed Tully-Fisher relation, this implies that the observed rotation 
speed of Tully-Fisher galaxies is within ~ 10 per cent of the typical circular speed of 
haloes hosting brightest galaxies of the same luminosity. 

Key words: galaxies: groups - galaxies: haloes - galaxies: clusters: general - large- 
scale structure of Universe. 



1 INTRODUCTION 

The abundance of galaxy clusters provides a very direct way 
to estimate erg, the linear theory rms density fluctuations in 
comoving spheres of 8/i -1 Mpc, extrapolated to redshift zero 
(e.g. Peebles, Daly & Juszkiewicz 1989; Frenk et al. 1990; 
White, Efstathiou & Frenk 1993; Eke, Cole & Frenk 1996; 
Viana & Liddle 1996; Ikebe et al. 2002; Pierpaoli et al. 2003; 
Schuecker et al. 2003; Henry 2004). This parameter sets the 
normalisation of the matter power spectrum. As the clus- 
ter abundance varies rapidly with as, the main difficulty 
in this method is not the determination of the abundance, 
but knowing the mass of the objects under consideration. 
The current ~ 10 per cent uncertainty on the value of ug 
(see Viana et al. 2003 and Henry 2004 for recent discussions 
of cluster-based estimates) brackets important differences 
in the formation histories of dark matter haloes. This has 
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implications for a broad range of topics, such as the red- 
shift at which gravitationally bound structures first form, 
the merger histories of galaxy-sized haloes and the concen- 
tration of the dark matter density profile within haloes. 

Many of the cluster samples used so far to determine 
ag have been based on the ROSAT X-ray All Sky Survey. 
It is clearly desirable to supplement these estimates with 
others derived from independent cluster samples, where the 
clusters are found using different techniques. There are al- 
ready a number of studies that have estimated the galaxy 
cluster mass function using catalogues of optically selected 
objects (e.g. Bahcall & Cen 1993; Biviano et al. 1993; Gi- 
rardi et al. 1998; Martinez et al. 2002; Bahcall et al. 2003). In 
this study, the 2dFGRS Percolation-Inferred Galaxy Group 
(2PIGG) catalogue (Eke et al. 2004a) is used, in conjunction 
with mock galaxy catalogues, to provide a new estimate of 
the mass function. The use of mock catalogues extends the 
previous work by quantifying some of the systematic biases 
inherent in the method. 
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The group luminosity function, that is the abundance 
of groups as a function of their total luminosity, has previ- 
ously been determined using galaxy redshift survey data by 
Moore, Frenk & White (1993), Marinoni, Hudson & Giuricin 
(2002) and Martinez et al. (2002). The 2PIGG catalogue in- 
creases the number of available groups by a factor of approx- 
imately two relative to the work of Martinez et al. (2002), 
who used the earlier data release from the 2dFGRS. The 
group luminosity function provides an indirect way to mea- 
sure the mass function if one knows how to relate mass to 
luminosity. For instance, if the statistical uncertainties on 
the cluster dynamical mass estimates are significantly larger 
than those on the cluster luminosities, then the combina- 
tion of the group luminosity function with a typical cluster 
mass-to-light ratio will yield a more accurate estimation of 
the cluster mass function than would be found from directly 
using the dynamical masses. Turning this around, given a 
theoretical mass function and a measured group luminosity 
function, one can infer the mapping from mass to luminosity 
as a function of group size (Marinoni & Hudson 2002). This 
picks out the group scales at which mass is most efficiently 
turned into starlight and, as such, encodes important clues 
about the process of galaxy formation (Benson et al. 2000; 
van den Bosch, Yang & Mo 2003). 

In this paper, mock galaxy catalogues, constructed from 
N-body simulations of the ACDM cosmology combined with 
a semi-analytical model for placing galaxies into dark matter 
haloes, are used to estimate how group-finding in the 2dF- 
GRS leads to biases in the recovered cluster mass and group 
luminosity functions. Then the 2PIGG catalogue is used to 
estimate the cluster mass and group luminosity functions, 
and the mass-to-light ratio for groups. This final quantity is 
inferred under the assumption that the group mass function 
is given by the fit to the halo mass function measured in N- 
body simulations, as described by Jenkins et al. (2001, J01). 
Given the typical halo mass, or equivalently circular speed, 
at a given luminosity, one can compare these results with 
the observed Tully-Fisher relation (Tully & Fisher 1977; 
Bell & de Jong 2001). Under the assumption that observed 
Tully-Fisher galaxies lie in typical haloes, this determines 
the conversion from observed rotation speeds of galaxies to 
halo circular velocities. This result has some bearing on the 
long-running debate concerning the ability of galaxy forma- 
tion models to match the galaxy luminosity function and the 
Tully-Fisher relation simultaneously (Kauffmann, White & 
Guiderdoni 1993; Cole et al. 1994; Heyl et al. 1995). 

Section 2 contains a brief description of the mock and 
real catalogues used in this study. The mass function calcula- 
tion and the constraint this imposes upon as is described in 
Section 3. Section 4 contains the group luminosity function 
results, which are applied to determining the group mass- 
to-light ratio variation in Section 5. In Section 6, the halo 
mass-to-light ratios are compared with the observed Tully- 
Fisher relation. 



2 BRIEF DESCRIPTION OF THE MOCK AND 
REAL CATALOGUES 

The construction of the 2PIGG catalogue is described in 
detail by Eke et al. (2004a). This paper also discusses the 
generation of mock galaxy catalogues from a combination of 
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Figure 1. The ratio of observed to total group luminosity for each 
mock group as a function of redshift. The total group luminosity 
is estimated using the global galaxy luminosity Schechter function 
described in the text to correct for galaxies outside the 2dFGRS 
flux limits. 

dark matter N-body simulations and semi-analytical galaxy 
formation models. The catalogue has been constructed from 
the two contiguous patches in the 2dFGRS (Colless et al. 
2001), which contain a total of ~ 190 000 galaxies. Of these, 
about 55 per cent are placed into ~ 29 000 groups containing 
at least two members. The median redshift of the groups 
is 0.11, like that of the galaxies, and the reliability of the 
estimated group masses and luminosities has been gauged 
using mock catalogues as described by Eke et al. (2004b). 

Mocks have been created from N-body simulations of 
the standard ACDM model with parameter values £7 m = 0.3, 
Q A = 0.7 and a 8 = 0.90 and 0.71. The a 8 = 0.90 simula- 
tion is the ACDM2 run described by Jenkins et al. (1998), 
whereas the lower as simulation cube contains 288 3 particles 
in a box of length 154/j" 1 Mpc, corresponding to a very simi- 
lar mass resolution. Both the semi-analytical models of Cole 
et al. (2000) and Benson et al. (2003) have been used to place 
galaxies into the dark matter distributions in the two simula- 
tion cubes. Note that the luminosities of the semi-analytical 
mock galaxies have been scaled by small, colour- preserving 
amounts so that the mock galaxy 6j-band luminosity func- 
tion exactly matches that found in the 2dFGRS. Dynamical 
masses of groups are estimated using 

M = A^, (2.1) 

where A = 5 (Eke et al. 2004a) , a is the 1-dimensional veloc- 
ity dispersion, calculated using the gapper algorithm (Beers, 
Flynn & Gebhardt 1990) and removing 85kms _1 in quadra- 
ture to account for redshift measurement errors, and r is the 
r.m.s. projected separation of galaxies from the group cen- 
tre, assuming an fi m = 0.3, S7a = 0.7 cosmological model. 
Observed group luminosities are estimated by summing the 
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luminosities of the individual galaxies with their associated 
weights to account for spectroscopic incompleteness: 

Lobs = T^WiLi. (2.2) 

This observed group luminosity is then corrected for galax- 
ies lying outside the survey flux limits assuming a global 
Schechter galaxy luminosity function 

0(L)d£ = 0»(A) a CX p(_i_) g, (2.3) 

with(M,-51og 10 /i,a) = (-19.725,-1.18) (Ekeetal. 2004b) 
and Mq = 5.33 in the fej band. Fig. 1 shows the ratio of 
L b s to the total group luminosity, Lj, T . Note that there are 
upper and lower flux limits imposed on the galaxies in the 
2dFGRS, hence the shape traced out by the galaxy groups in 
Fig. 1, with the upper flux limit rejecting bright galaxies in 
the nearby groups, and the lower flux limit being important 
at higher redshifts. The stripes traced out by the points are 
a result of the distinct flux limits in different patches of the 
survey. 

In all of what follows, the group sample will be restricted 
to the z < 0.12 objects, because at higher redshifts, less 
than half of the total group luminosity is contained in visible 
galaxies and the correction for this missing light becomes too 
uncertain (see Fig. 1). For the same reason, the volume at 
z < 0.02 will also be excised. 

When calculating the abundance of groups, it is neces- 
sary to know the volume surveyed. The maximum volume 
within which a particular group could have been detected by 
a survey is Knax ■ 

If a group is defined to comprise at least N 
galaxies, then the group Knax could be approximated by the 
Knax of the JVth most luminous galaxy in that group. This 
takes into account the variable sky coverage as a function 
of flux in the 2dFGRS. However, for N > 1, this ignores 
the possibility that the set of galaxies linked together by the 
friends-of-friends (FOF) algorithm may not remain grouped 
at different redshifts, when the linking volume has changed 
in size and the set of galaxies satisfying the flux limits of 
the survey may have changed to the extent that parts of 
the group are no longer joined. To deal with these issues 
exactly would be awkward. These complications are minor 
though, because the linking volume changes with redshift to 
take into account the varying sampling of the galaxy distri- 
bution. Consequently, the Knax of the Nth most luminous 
galaxy is used to define the group Knax throughout the rest 
of this paper. The mock catalogues can be used to gauge 
the extent to which this choice biases the results. Of course, 
for the group luminosity function, iV m i n = 1 can be used, in 
which case the K max value is correct. 



3 CLUSTER MASS FUNCTIONS 

The accuracy with which dynamical masses can be esti- 
mated varies rapidly with the number of galaxies sampling 
the group potential (see figure 3 of Eke et al. 2004b). Er- 
rors in estimated masses coupled with the steep decline in 
the abundance of clusters, i.e. groups with masses above 
~ 10 14 /i — 1 Mq, as a function of increasing mass, lead to sys- 
tematic biases in the inferred abundance of clusters. This 
happens because more low-mass objects are scattered to 
higher masses than are scattered down from those higher 
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Figure 2. Model and mock recovered mass functions. The solid 
line traces the mass function of dark matter haloes in the simu- 
lation cube (the ACDM2 run by Jenkins ct al. 1998). The other 
lines show the mass function inferred using different minimum 
numbers of galaxies to define the group sample, as indicated in 
the key. Statistical errors, shown by the error bars, are calculated 
using the scatter between 22 mock catalogues constructed from 
the Hubble Volume simulation. 



masses. Consequently, the abundance of high mass galaxy 
clusters will be systematically overestimated as a result of 
the uncertainties in measuring the individual cluster masses. 
This effect is shown in Fig. 2 for the mock catalogues, which 
compares the mass functions recovered from groups samples 
of different AT m i n at z < 0.12 in the mocks with the mass 
function of dark matter haloes in the simulation cube. The 
model mass function shows the abundance of dark matter 
haloes found using a friends-of-friends (FOF) group-finder 
with a linking length b — 0.2 times the mean interparticle 
separation (Davis et al. 1985). Error bars represent the stan- 
dard deviation among the results from 22 different mock cat- 
alogues made using the Hubble Volume simulation (Evrard 
et al. 2002). Norberg et al. (2002) provide a detailed descrip- 
tion of how the mock catalogues are constructed. Hawkins 
et al. (2003) have shown that the galaxies in mock cat- 
alogues made from the Hubble Volume simulation cluster 
like those in the 2dFGRS, suggesting that the use of these 
mocks is appropriate. Cosmic variance only increases the 
statistical errors by ~ 10 per cent relative to the 1/Knax- 
weighted Poisson errors. Note that restricting the sample 
to include only groups containing large numbers of galaxies 
introduces incompleteness at lower masses, which, to some 
extent, counteracts the overestimation of the mass function 
due to mass measurement errors. Nevertheless, as Fig. 2 
shows, the abundance is still overestimated by a factor of 
~ 2 for M > 3 x 10 14 /i _1 Mq. For the remainder of this Sec- 
tion, z max = 0.12 and iV m i n = 20 will be used. These cuts 
restrict the sample to ~ 350 clusters. 

The dynamically-inferred group masses show more scat- 
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Figure 3. Mass functions for dark haloes in the simulation and 
for groups identified in the mock catalogue. A solid line traces 
the mass function of dark matter haloes in the simulation cube. 
Also plotted are mass functions for mock groups inferred using 
dynamical masses (short-dashed line) and group luminosities with 
a shift for the typical mass-to-light ratio (dotted line). The long- 
dashed line represents the J01 fitting formula to describe the dark 
matter halo mass function appropriate for the power spectrum of 
the model. 



ter about the true halo mass in the simulations than masses 
inferred using the total group luminosity and a typical mass- 
to- light ratio (Eke et al. 2004b). Thus, one can produce a 
more robust estimate of the cluster mass function by using 
the cluster total luminosity function and a typical mass-to- 
light ratio. This is shown, at least for the highest masses, in 
Fig. 3. The median mass-to-light ratio for all mock groups 
with L > W ll r 'h~ 2 L Q (T = 471 hM Q /L Q ) has been used as 
a global shift to the cluster luminosity function to infer the 
mass function (dotted line). The mass-to- light ratio varies 
little for these most massive objects (Eke et al. 2004b), so 
this simple shift is appropriate. Also shown is the J01 fitting 
formula for the mass function that describes the population 
of haloes from which the model curve is sampled. Hereafter 
in this paper, only recovered mass functions obtained using 
the cluster luminosity function and the cluster mass-to- light 
ratio, as shown by the dotted line in Fig. 3, will be consid- 
ered. 

In order to place constraints upon the normalisation of 
the power spectrum, it is necessary firstly to model the bias 
in recovering the mass function and secondly to compare 
the measured mass function in the 2PIGG catalogue with 
the biased models for each value of erg. Fig. 4 shows the 
J01 mass functions for two different values of as, before and 
after applying a shift to account for the measurement errors. 
An approximation to the bias introduced by measuring the 
mass function is to translate the J01 curves to the right by 
0.17 and 0.20 for ag = 0.90 and 0.71 respectively. These 
shifts provide the best \ 2 fits (assuming weighted Poisson 



Figure 4. The erg dependence of the recovered mass function. 
Long-dashed lines show the J01 fitting function describing the 
mass function for the two different values of erg used to make mock 
catalogues. The solid lines show the effect of measurement errors 
on these curves, as described in the text. A dotted line traces the 
mass function recovered from the erg = 0.9 mock, whereas the 
short-dashed line is the corresponding curve for erg = 0.71. The 
horizontal arrows show the impact of the measurement bias. 



errors) to the mass functions measured from the mock group 
catalogues for both values of ag. The uncertainty on the size 
of these translations, judging by the change in \ 2 of the 
resulting fits, is ~ 7 per cent in mass for the Cole et al. 
(2000) mocks. The shifts required for the mocks made using 
the Benson et al. (2003) semi-analytical procedure differ by 
~ 7 per cent from the Cole et al. (2000) mocks. Thus, one 
should consider these corrections to be uncertain by ~ 10 
per cent in mass. With these shifts, each model curve can be 
converted into a mass function in measurement space. Note 
how the measurement bias is larger for the lower value of 
ag, where the clusters are more difficult to detect because of 
the lower contrast. This, unfortunately, reduces the power 
of the test, but it is nevertheless apparent that these two 
different values of as can still be distinguished. 

The mass functions inferred from the 2PIGG data are 
shown in Fig. 5. Once again, the directly determined mass 
function, based on dynamical masses, gives higher abun- 
dances at this mass than the mass function recovered from 
the luminosity function and the typical cluster mass-to-light 
ratio (T = 429 hM e /L Q ). Note that the 2PIGG results gen- 
erally lie between the convolved model curves for as = 0.71 
and 0.90 (fi m = 0.3 in both cases). Error bars again are the 
1/lmax-weighted Poisson errors. This neglects the ~ 10 per 
cent additional contribution from cosmic variance, which is 
minor relative to the systematic uncertainties inherent in 
this method. 

By assuming that the convolving function generates a 
mean shift that varies linearly with as and not at all with 
f2 m , one can extrapolate the convolved model to a range 
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Figure 5. 2PIGG and model group mass functions. The model 
mass functions after convolution with the measurement errors, 
shown with solid lines in the previous figure, are now shown with 
dotted (erg = 0.90) and dashed (erg = 0.71) lines. The points with 
error bars show the 2PIGG results using the cluster luminosity 
function and typical mass-to-light ratio. A solid line shows the 
2PIGG mass function based on the dynamical mass measurements 
of equation 2.1. 
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Figure 6. The constraints on as and f2 m from the 2PIGG cata- 
logue. The x 2 contours are shown, with (solid lines) and without 
(dotted lines) the correction for the bias in the recovery of the 
mass function. The contours contain 68, 95 and 99.7 per cent 
of probability, and a cross marks the best-fitting parameter set 
for each case. The dashed line shows a fit to the most probable 
o"8(f2m) for the bias-corrected contours. 



Figure 7. Comparison of the 2PIGG results (solid black line) 
with various other estimates of the amplitude of matter fluctu- 
ations. The red lines show measurements based largely on the 
cluster abundance: Eke, Cole & Frenk (1996, solid line); Pier- 
paoli et al. (2001, short dashed); Bahcall et al. (2003, short-long 
dashed); Reiprich & Bohringer (2002, dot-short dashed); Seljak 
(2002, dotted); Viana, Nichol & Liddle (2002, long dashed); Allen 
et al. (2003, dot-long dashed). Blue lines trace the amplitudes 
inferred using weak lensing-based estimates from: Brown et al. 
(2003, dotted); Rhodes et al. (2004, long dashed); Heymans et al. 
(2005, dot-short dashed); Massey et al. (2005, solid); van Waer- 
beke, Mellier & Hoekstra (2005, short dashed). The green line 
shows the results of Tinker et al. (in preparation), who use a 
model of the halo occupation distribution, in conjunction with 
the galaxy correlation functions for different luminosity galaxies 
in the 2dFGRS and the 2PIGG mass-to-light ratio variation with 
group size. The rounded 68 and 95 per cent probability contours 
are from figure 12 of Sanchez et al. (2005), and result from a joint 
analysis of the power spectra of both the CMB and 2dFGRS. 

of $l m and as, in order to determine the best-fitting val- 
ues of these parameters. The results of this exercise are 
shown in Fig. 6. Only the 5 data points with 14.4 < 
log[Af/(/i _1 M Q )] < 15 are used in the \ 2 nt - Tne contours 
contain 68, 95 and 99.7 per cent of the probability. Dotted 
and solid lines show the constraints on the parameters before 
and after the correction for the measurement bias is applied. 
In short, this correction changes the best-fitting parameters 
from (fi m = 0.39, as = 0.88) to (0.25, 0.90). The dashed line 
showing the variation of the best-fitting ag as a function of 
S7 m in the range 0.18 < f2 m < 0.50 has the equation: 

a 8 = 0.25fi-°- 92 - 4 ' 5(n ™- a22)2 . (3.1) 

Note that these contours have all been calculated under the 
assumption that the CDM transfer function has a shape 
given by F — 0.15 (Bardeen et al. 1986; Sugiyama 1995) as 
has recently been measured using the 2dFGRS (Cole et al. 
2005). Using values of T = 0.1 or 0.25 changes the value of 
as by less than 5 per cent for f2 m = 0.3. While the statistical 
uncertainty upon the estimate of Cg(f2 m ) is small (~ 10 per 
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cent), the systematic correction that has been applied to ac- 
count for the measurement bias is substantially bigger (~ 20 
per cent). It is worthwhile reiterating that the probability 
contours involve extrapolating the noisy fits of convolving 
functions to the mock recovered mass functions, so this sys- 
tematic shift is quite rough. As was discussed above, the 
correction is uncertain by ~ 10 per cent in mass. This corre- 
sponds roughly to 5 per cent systematic uncertainty in the 
estimated as value from this shift. There is another system- 
atic uncertainty associated with the calibration of the mass 
estimation (the value of A in equation 2.1). 

An illustration of the relative importance of system- 
atic uncertainties in current estimates of as is contained in 
Fig. 7, which shows the results from a number of different 
studies. The 2PIGG line is shown in black, along with other 
cluster abundance-based estimates in red, and weak lensing- 
based results in blue. The black probability contours show 
the results of Sanchez et al. (2005) , who jointly analysed the 
CMB and 2dFGRS power spectra. The green line shows the 
best-fitting parameters inferred using the 2dFGRS correla- 
tion functions, a model for the halo occupation distribution 
and the 2PIGG mass-to- light ratios (Tinker et al., private 
communication). Quoted errors on the various curves are 
usually of the order of 10 — 15 per cent. Some of the vari- 
ation from one study to another will come from assuming 
different power spectrum shapes (F), but this dependence 
is quite weak (Seljak 2002; Rhodes et al. 2004) so the ex- 
cess scatter between the curves is a sign of other systematic 
effects. For the cluster abundance-based measurements, the 
main source of systematic differences comes due to the choice 
of the mass-observable relation (see Henry 2004). 



4 GROUP TOTAL LUMINOSITY FUNCTION 

The total &j-band group luminosity function is a distribution 
that, unlike the galaxy luminosity function, can readily be 
compared with theoretical dark matter halo mass functions. 
In Section 5 the group luminosity function will be used in 
conjunction with theoretical mass functions to derive halo 
mass-to-light ratios. 

Fig. 8 shows how well the 6j-band group luminosity 
function can be recovered from a mock catalogue. Only the 
groups with 0.02 < z < 0.07 are included in the calculation 
of these curves. As the upper limit in redshift is increased, 
the typical correction for luminosity in galaxies fainter than 
the flux limit increases. Consequently, the recovered lumi- 
nosity functions become increasingly biased high. On the 
other hand, if the maximum redshift is too restrictive, then 
there are insufficient clusters, and the volume probed may 
no longer be representative of the Universe as a whole. The 
choice of z max = 0.07 is a compromise between these two 
competing factors. Even with this value, the group luminos- 
ity function recovered from the mock catalogue is still biased 
slightly high relative to the true model function for the most 
luminous systems in the simulation. Note how the value of 
iV m in impacts upon the abundance of low luminosity groups. 
The best recovery of the model occurs when iV m i n = 1. Some 
depletion of the less luminous groups is evident: these ob- 
jects sometimes contaminate the bigger groups, causing the 
abundance of the more luminous systems to be slightly over- 
estimated. 
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Figure 8. The abundance of groups as a function of their total 
fej-band luminosity. A solid line traces the actual group luminosity 
function in the simulation, and the other lines show the recovered 
functions using all mock groups with z < 0.07 and N > iV m i n . 
Error bars are shown only for the N mln = 10 case, and are cal- 
culated using the scatter between 22 different mocks made from 
the Hubble Volume simulation. 



Fig. 9 shows how the group luminosity function recov- 
ered from the 2PIGG catalogue compares with that in a 
mock catalogue constructed from a simulation with as = 0.9 
using the Cole et al. (2000) semi-analytical model. Only 
groups at 2 < 0.07 that contain at least one detectable 
galaxy are included in the calculation. The 2PIGG results 
accurately trace those recovered from the mock. Also shown 
is the 2dFGRS galaxy luminosity function, which crosses 
the group luminosity function at ~ L», the characteristic 
luminosity in the galaxy luminosity function (Norberg et al. 
2002). Many of the fainter galaxies reside in groups with 
L > L* ; hence the abundance of low luminosity groups lies 
beneath that of the galaxies. 

Fig. 10 shows how the 2PIGG results compare with 
other published work. The group luminosity functions of 
Marinoni et al. (2002), from the Nearby Optical Galaxy 
(NOG) catalogue, Moore, Frenk & White (1993), from the 
CfA survey, and Martinez et al. (2002), using the 2dF 
Galaxy Group Catalogue (2dFGGC), are all quite simi- 
lar to that from the 2PIGGs for group luminosities above 
~ 10 n /i~ 2 L Q . At lower luminosities, the 2dFGGC results 
lie lower, because of the requirement that their groups con- 
tain at least 4 galaxies. The 2PIGG group luminosity func- 
tion also lies somewhat above those of Moore et al. and 
Marinoni et al. Note that the Marinoni et al. curve has 
been shifted 0.55 magnitudes fainter from their published 
fit to account for both the 0.3 magnitudes difference be- 
tween their B band (RC3 asymptotic photometric system) 
and 6j (Marinoni, private communication), and 0.25 mag- 
nitudes of internal absorption that they had included. The 
Moore et al. curve has been shifted under the assumption 
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Figure 9. The abundance of groups as a function of their total 
bj-band luminosity. A solid line traces the model group luminosity 
function, the dot-dashed line shows the recovered function using 
all mock groups with z < 0.07 and N > 1. The points show the 
group luminosity function found using the 2PIGG catalogue. 
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Figure 10. The abundance of groups as a function of their total 
fej-band luminosity. The points show the group luminosity func- 
tion found using the 2PIGG catalogue out to z = 0.07, whereas 
the lines show the results of Moore, Frenk & White (1993, solid), 
Marinoni, Hudson & Giuricin (2002, dashed) and Martinez et al. 
(2002, dotted). 



that &j = -Bzwicky — 0.05. This should make these curves di- 
rectly comparable with the other results in the bj band. Both 
the Moore et al. and Marinoni et al. studies probe smaller 
volumes than in the 2PIGG case, so cosmic variance may be 
partly responsible for the lower abundance of <J 10 11 hT 2 L Q 
groups found in these earlier studies. Also, the group-finder 
used by Marinoni et al. places a higher fraction of galaxies 
into groups, so one might expect that fewer isolated galaxies 
would remain, yielding fewer low luminosity groups. 



5 GROUP MASS-TO-LIGHT RATIOS 

By comparing, at a fixed abundance, the measured cumula- 
tive group luminosity function with a mass function moti- 
vated by numerical simulations of a CDM model, it is pos- 
sible to determine the mass-to-light ratio down to groups 
that contain only a single visible galaxy in the 2dFGRS. 
This method (see Marinoni & Hudson 2002, MH02) is shown 
schematically in Fig. 11, where the mass-to- light ratio gives 
the mapping from the cumulative luminosity function to the 
cumulative mass function. An assumption of this method is 
that the group mass varies monotonically with luminosity, so 
that this mapping is unique. Fig. 11 shows the importance 
of the choice of mass function. For instance, the Press & 
Schechter (1974) formula underestimates the abundance of 
the largest clusters and overestimates that of Local Group- 
sized objects. This overabundance would produce mass-to- 
light ratios for small groups that are too high by ~ 50 per 
cent. 

This is an indirect, model-dependent measurement of 
the mass-to- light ratio, because the mass function shape and 
amplitude are assumed to be of a particular form, rather 
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Figure 11. Cumulative group abundances as functions of both 
mass and luminosity in the simulation model. The quantity X 
on the horizontal axis, is cither the group bj-band luminosity in 
h~ 2 Lq or the group mass in ft. -1 Mq. A dotted line shows the 
J01 fit to the mass function, whereas the dashed line is the corre- 
sponding Prcss-Schcchter curve. The horizontal distance with an 
arrow represents the inferred group mass-to-light ratio for groups 
with a luminosity equal to 10 11 /i _2 Lq. 
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Figure 12. Model mass-to-light ratios as functions of group lu- 
minosity. The solid curve traces the median mass-to-light ratio in 
bins of group luminosity, whereas the long dashed line represents 
the function required to map the model group luminosity func- 
tion to the model group mass function. If, instead of using the 
actual group mass function from the simulation, the J01 or Press- 
Schcchtcr formulae are used, then the dotted and short dashed 
curves result. 
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Figure 13. The recovery of the mass-to-light ratio as a function 
of group luminosity. The dotted curve traces the median mass-to- 
light ratio recovered from the mock catalogue using the dynam- 
ically inferred masses. This is the attempt to recover the model 
represented by the solid curve. The three sets of points show the 
mass-to-light ratio inferred from different sets of mock groups 
using the abundance-matching technique, and arc attempts to 
measure the long-dashed curve of the model. 



than being directly measured. However, with this method 
one can probe to smaller group sizes, where group lumi- 
nosities are still well-defined but the masses are difficult to 
measure. The uncertainty associated with the accuracy of 
the adopted mass function includes a non-negligible con- 
tribution from uncertainty in the value of as for the real 
Universe. As discussed in Section 3, the abundance of such 
clusters depends very sensitively upon as- Consequently, the 
mass-to- light ratios inferred from this method for these large 
objects will suffer from this uncertainty. However, for smaller 
groups, this effect is minor, and the mass function is quite 
robust to changes in as- Thus, this method for measuring 
the mass-to-light ratio is complementary to the direct mea- 
surement technique, which is most effective for the larger 
systems that have enough galaxies to allow an accurate mass 
estimation. 

It is worth noting that the mass-to- light ratio as a func- 
tion of group luminosity obtained by matching abundances 
is not the same as the median mass-to-light ratio of groups 
in each particular luminosity bin. This is because there is 
some scatter in the group masses at each group luminosity. 
However, the difference between these two mass-to-light ra- 
tios in the model can be seen in Fig. 12 to be small. Also 
shown are the results of inferring the mass-to-light ratio as- 
suming the J01 and Press-Schechter mass functions rather 
than using the actual mass function. This illustrates the 50 
per cent overestimation for small haloes mentioned in the 
discussion of the preceding figure. The J01 mass function is 
used throughout the rest of this paper. 



5.1 Results 

The systematic errors in the recovery of the group lumi- 
nosity function alluded to in Section 4 impact significantly 
upon the inferred variation of the group mass-to-light ratio. 
Fig. 13 shows how changing the maximum redshift of the 
mock group sample used to measure the luminosity func- 
tion affects the results. For the largest groups, the mass-to- 
light ratio drops with increasing z max , as the abundance of 
clusters becomes more overestimated because of increasing 
contamination by interlopers and larger corrections for lu- 
minosity in galaxies beneath the flux limit. For groups with 
L <. 10 10 /i~ 2 L however, the volume probed is smaller and 
the errors associated with the correction from observed to to- 
tal group luminosity are smaller. Consequently, the group lu- 
minosity functions are less dependent upon z m ax- The direct 
method for inferring the group mass-to-light ratio becomes 
significantly biased for L < 10 10 /i -2 Lq, so this group lumi- 
nosity is an appropriate value at which to switch between 
the direct and indirect methods. In this way, the 'corrected' 
mass-to-light ratio can be recovered over almost four orders 
of magnitude in group luminosity. For L > 10 10 h~ 2 L the 
appropriate correction factor is the difference between the 
directly recovered curve (dotted) in Fig. 13 and the model 
curve (solid). For L < 10 10 /i~ 2 Lq the correction factor is 
chosen to be the difference between the crosses and the 
model shown by the solid curve, so that the final 'corrected' 
curve provides an estimate of the typical mass-to- light ratio 
as a function of group luminosity. 

The uncorrected results for the 2PIGG groups are 
shown in Fig. 14 and compared with those of MH02. As for 
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Figure 14. The recovery of the mass-to-light ratio as a func- 
tion of group luminosity for the 2PIGG sample. The dotted curve 
traces the median mass-to-light ratio recovered using the dynam- 
ically inferred masses. The three sets of points show the mass-to- 
light ratio variation inferred from different sets of 2PIGG groups 
using the abundance-matching technique. A dashed line shows 
the results of Marinoni & Hudson (2002) from the NOG sample 
of galaxy groups. 



Figure 15. The 'corrected' mass-to-light ratio as a function of 
group luminosity for the 2PIGG sample. The error bars include 
the statistical errors added in quadrature to the uncertainty in the 
systematic shift applied to 'correct' the measurements. Directly 
measured mass-to-light ratios are shown with points, whereas 
those inferred assuming that the global halo mass function is given 
by the fitting function of J01 are shown by the shaded region. 



the mock catalogues, the systematic differences between the 
different z max samples are evident for the bigger groups. The 
increasing overestimation of the group luminosity function, 
as z max increases, gives rise to decreasing inferred mass-to- 
light ratios. Note that the inferred mass-to-light ratio for the 
smallest groups is comparable with that measured directly 
for the clusters. 

It is apparent that the 2PIGG results are somewhat 
different from those of MH02. This discrepancy arises for 
two main reasons. For the small groups, MH02 use the 
Press-Schechter formula for the mass function. As was seen 
in Fig. 13, this leads to an overestimation of the mass- 
to-light ratio of small groups by ~ 25 per cent (this, of 
course, assumes that the underlying mass function is that 
appropriate for ACDM). This difference is compounded at 
L bj < 10 10 - 6 /i~ 2 L Q by the different group luminosity func- 
tions measured from the 2PIGG and NOG samples, as was 
shown in Fig. 10. As the NOG sample contains a lower abun- 
dance of groups at L < 10 10 6 /i _2 L Q , the inferred mass-to- 
light ratio is larger than for the 2PIGG sample. Note that 
the MH02 results in Fig. 14 differ from those plotted in fig- 
ure 15 by Eke et al. (2004b), because they did not use the 
appropriate waveband correction. 

Having calibrated the biases introduced by the two 
methods to measure group mass-to-light ratios with the 
mock catalogues, it is now possible to 'correct' the 2PIGG 
results in Fig. 14. Using the 2 max = 0.07 data for the indi- 
rect mass-to- light ratio measurement below L = 10 10 /i~ 2 L 
yields the results shown in Fig. 15. The errors on the indi- 
rectly measured points come from the scatter between 22 



mock catalogues made from the Hubble Volume simulation, 
added in quadrature to the uncertainty in the correction fac- 
tor. For the L > 10 10 /i~ 2 L results, the error bars show the 
statistical uncertainty in the raw mass-to-light ratio added 
in quadrature to the uncertainty in the correction factor. 
The uncertainties in the correction factors are calculated 
by taking eight different mock catalogues with different as 
values and semi-analytical schemes, and finding the scatter 
between the corrections. 

The corrected mass-to-light ratio variation with lumi- 
nosity for the 2PIGGs shows a plateau at cluster masses, 
a decrease to a minimum at L ~ 10 10 /t~ 2 LQ, and an in- 
crease for smaller groups. The value of the halo mass at 
this minimum is ~ lO 12 ^ 1 Mq, which agrees well with the 
weak lensing analysis of Hoekstra, Yee & Gladders (2004). 
They found a halo of this luminosity to have M200 = 
(8.4 ± 0.7 ± 0.4) x lO 11 /^ 1 M Q , which should be increased 
by ~ 20 per cent to compare with the halo mass definition 
used here. It is reassuring that these two completely different 
methods show such consistency. 

One might wonder how surprising it is that the cor- 
rected, recovered mass-to-light variation looks so similar to 
that in the model for L < 10 10 /i~ 2 Lq, where the masses 
of the 2PIGGs are not directly measured. After all, the 
model mass function has been assumed to be appropriate 
for the real Universe, and the galaxy luminosity function in 
the model has been scaled so as to match that of the 2dF- 
GRS. The only possibility for the corrected 2PIGG results 
to differ from the model in Fig. 15 is if the group lumi- 
nosity functions differ, despite the identical global galaxy 
luminosity functions. Thus, the measurement of an increase 
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Figure 16. The ratio of the median brightest galaxy luminos- 
ity to the total group luminosity as a function of luminosity for 
the mock groups and 2PIGGs, shown with solid and dotted lines 
respectively. 



in 2PIGG mass-to-light ratio as group size decreases from 
L ~ lO lo /i~ 2 L/0 is, if not inevitable, to be expected merely 
from the fact that the faint end slope of the galaxy luminos- 
ity function is flatter than that of the ACDM mass function 
at low masses. The extra information available here comes 
from the fact that the group luminosity function is tracing 
the same structures probed with the mass function measure- 
ments. Consequently, the normalisation of the mass-to-light 
ratio is not predetermined but, in fact, recovered in addition 
to the variation with group size. 



6 THE TULLY-FISHER RELATION 

The results in Fig. 15 show how the typical halo mass varies 
with halo luminosity. For a particular definition of a halo, 
one can convert this to a relation between the halo circular 
speed and the halo luminosity via 



L b , = 



47rAp c 



(6.1) 



where Ap cr i t represents the mean enclosed density of the 
halo, taken to be 100 times the critical value. (Note that this 
value is not exact because the identification of the groups is 
done using a FOF algorithm, rather than growing a sphere 
out to a particular density contrast.) This relation is reminis- 
cent of the Tully-Fisher relation between galaxy luminosity 
and rotation speed at the edge of the visible galaxy. For low 
luminosity systems, a central bright galaxy usually domi- 
nates the total group luminosity, so the galaxy and group 
luminosities are likely to be quite similar. Fig. 16 shows the 
median ratio of the brightest galaxy to total group luminos- 
ity, as a function of total group luminosity. The mock and 
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Figure 17. The galaxy rotation speed, Ufl at , and inferred halo 
circular speed, v c , as functions of blue luminosity. For the galaxy 
rotation speed curve (dotted line), the luminosity is that of the 
galaxy, according to the Tully-Fisher relation of Bell and de Jong 
(2001). The luminosity for the dashed curve is the total group 
luminosity multiplied by the ratio shown in Fig. 16 to convert it 
to the typical luminosity of the brightest galaxy in each halo. The 
solid line shows the ratio of these two velocities as a function of 
primary galaxy luminosity. 



2PIGG curves are very similar, with the main difference be- 
ing that the brightest galaxy in the bigger mock groups typ- 
ically contains a slightly higher fraction of the total group 
luminosity than in the 2PIGGs. 

If one were to apply the median shift of Fig. 16 to con- 
vert the group luminosity to the typical luminosity of the 
brightest contained galaxy, then equation (6.1) could relate 
galaxy luminosity to halo circular speed at the virial radius. 
Then, one might reasonably ask, given the Tully-Fisher re- 
lation and the results in Fig. 15, what is the relationship 
between the galaxy rotation speed and the halo circular 
speed for galaxies with L < 10 10 /i~ 2 Lq? The answer to 
this question is contained in Fig. 17. The B-band Tully- 
Fisher relation of Bell & de Jong (2001) has been converted 
to 6j, and is shown with a dotted line. Note that these au- 
thors advocate the use of vaat to characterise the velocity 
derived from the galaxy rotation curve. The dashed curve 
shows the halo circular velocity from equation (6.1) as a 
function of luminosity (ie the group luminosity multiplied 
by the factor shown in Fig. 16 to convert it to a galaxy lu- 
minosity). A solid line traces the ratio of ug at to v c . How 
should this figure be interpreted? If one considers a galaxy 
with L 6j ~ 10 10 /i~ 2 L Q , then the Tully-Fisher relation of 
Bell & de Jong gives «fl at ~ 160 km s -1 , whereas the abun- 
dance matching method implies that it typically lives in a 
halo with circular speed ~ 200 km s -1 . For less luminous 
galaxies, these two velocities become more similar, so that 
the flat part of the galaxy rotation curve corresponds to a 
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speed that is approximately the same as the circular velocity 
of the host halo. 

The suggestion that the galaxy rotation speed is sim- 
ilar to the halo circular velocity is broadly in accord with 
the results for the /-band Tully-Fisher relation in the semi- 
analytical galaxy formation models of Cole et al. (2000) and 
Croton et al. 2005. They found that the galaxy luminosity 
function and the Tully-Fisher relation could be matched si- 
multaneously when the galaxy rotation speed was assumed 
to equal the halo circular velocity. However, one would ex- 
pect as baryons concentrate in the centres of galaxy-sized 
potential wells, they would drag in dark matter, leading to 
galaxy rotation speeds that are a few tens of per cent larger 
than the halo circular velocity. In this case, the Tully-Fisher 
relation is no longer reproduced by the model. Similar con- 
clusions have been drawn from hydrodynamical simulations 
of galaxy formation by Navarro & Steinmetz (2000) and Eke, 
Navarro & Steinmetz (2001), and also from the analysis of 
de Jong et al. (2004). 

It is possible that the comparison above could be mis- 
leading; perhaps the galaxies included in Tully-Fisher sam- 
ples are atypical. If the observed set of galaxies lie in either 
unusually luminous or unusually low concentration haloes, 
then this might reconcile the apparent discrepancy between 
the low observed rotation speeds and the higher expected 
values. For example, the semi-analytical model behind the 
mock catalogues used in this work predicts that the distri- 
bution of masses of haloes with L ~ 10 10 /i~ 2 Lq has a la 
scatter that equates to ~ 20 per cent in halo circular speed. 
If, for some reason, the galaxies selected in the Tully-Fisher 
measurements lie in the less massive haloes at that particu- 
lar luminosity, then this could have a very significant impact 
in reducing the dark matter contribution to the observed ro- 
tation curves. Similarly, if typical Tully-Fisher galaxies hap- 
pen to lie in atypically unconcentrated haloes, then the dark 
matter contribution to the observed rotation curves would 
again be reduced. 



7 CONCLUSIONS 

The 2PIGG catalogue has been used to measure the mass 
and luminosity functions of groups and clusters. By com- 
bining the measured abundance of clusters as a function 
of luminosity, with a typical mass-to-light ratio, the clus- 
ter mass function is measured more accurately than by us- 
ing the dynamically-inferred masses directly. After removing 
a bias due to measurement errors in the abundance, using 
that found in mock catalogues constructed from ACDM N- 
body simulations, the 2PIGG mass function implies that 
erg = 0.25 fim - 92 - 4 - 5 («m-o.22) 2 for lg < flm < o.50. While 
the statistical uncertainty on this value is around 10 per 
cent, the systematic correction that has been applied is 
closer to 20 per cent. The uncertainty in this correction, 
coupled with the systematic uncertainty in the normalisa- 
tion of the estimated cluster masses (the value of A in equa- 
tion 2.1), currently limits the precision of this determination 
of ag to ~ 20 per cent. To decrease this uncertainty signifi- 
cantly would require a more detailed understanding of how 
the observed galaxies populate the underlying dark matter 
haloes. 

The group luminosity function is in good agreement 



with previous studies for groups with total fej-band lumi- 
nosity exceeding ~ 4 x 10 10 hT 2 Lq. At lower luminosities, 
the abundance of 2PIGGs is somewhat larger than has previ- 
ously been found from smaller samples. Matching this abun- 
dance to that of haloes of a given mass in ACDM simulations 
gives the mass-to-light ratio as a function of group size. This 
indirect method combined with the direct estimates of Eke 
et al. (2004b) at higher luminosities allows the recovery of 
the variation of the group mass-to-light ratio over almost 
four orders of magnitude in group luminosity. The resulting 
function has a minimum value of ~ 100 HMq/Lq at group 
luminosity L bj ~ 5 x 10 9 /i _2 Lq. Finally, the group mass- 
to-light ratio was used to infer the halo circular speed as a 
function of group luminosity. Comparison to the observed 
Tully-Fisher relation for galaxies suggests that, at a given 
galaxy luminosity, the observed rotation speed is similar to 
the typical halo circular speed. 

The dependence of the group mass-to-light ratio on 
group luminosity provides a key test of various feedback pro- 
cesses invoked in models of galaxy formation (Benson et al. 
2000; Benson et al. 2003; Croton et al. 2005). The agreement 
between the indirect determination at low luminosities pre- 
sented in this paper and the results of the semi-analytical 
ACDM model assumed in the construction of the mock cat- 
alogues is encouraging. A direct determination based, for 
example, on a deeper redshift survey is desirable. 
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